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Bs Mesons using Staggered Light Quarks 

Matthew Wingate,^* Junko Shigemitsu,^ G. Peter Lepage,*^ Christine Davies,'^ and Howard Trottier'^ 
"^Department of Physics, The Ohio State University, Columbus, Ohio 43210, USA 
''Newman Laboratory of Nuclear Studies, Cornell University, Ithaca, NY 14853, USA 
■^Department of Physics & Astronomy, University of Glasgow, Glasgow, G12 8QQ, UK 
'^Physics Department, Simon Fraser University, Burnaby, B.C., V5A 1S6, Canada 

Last year we proposed using staggered fermions as the light quarks, combined with nonrelativistic heavy quarks, 
in simulations of heavy-light mesons. A first round of tests which focuses on the Bs meson has been completed 
using quenched lattices, and results are presented here for the kinetic Bs mass, the B* — Bs splitting, and /b, ■ 
The next project, already underway, is to compute the B and Bs decay constants and spectra on the n/ =2+1 
and 3 MILC lattices. We report on progress with one set of these configurations. 



1. INTRODUCTION 

Although a variety of lattice methods have 
been used to simulate the heavy quark within a 
heavy-light meson, the light quark description has 
not been studied as thoroughly; only Wilson-like 
fermions have been used in large scale simulations 
to date. Recent simulations of light hadrons sug- 
gest, however, that improved staggered fermions 
are superior to Wilson-like fermions, particularly 
for small quark masses. Improved staggered 
fermions have an exact chiral symmetry at zero 
mass, no 0{a) errors and excellent scaling, and 
they require significantly less computational ex- 
pense. These benefits will facilitate the study 
of leptonic decay constants and form factors in 
semileptonic decays as the physical light quark 
limit is approached. 

We have been studying the coupling of stag- 
gered light quarks to nonrelativistic heavy quarks 
for simulation of heavy-light mesons. Since our 
initial study last year [1] we have resolved a ques- 
tion regarding the improved staggered action on 
a very coarse S'^ x 20 lattice and have completed 
a calculation of the Bg mass, B* — Bs mass split- 
ting, and on a finer 12'^ x 32 lattice. We 
present those results in Section 2. In Section 3 
we discuss progress with this simulation method 
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on the 2-1-1 flavor MILC configurations. 

Given the constraints of this write-up we refer 
the reader to [1] for a reminder of our method 
for constructing operators out of NRQCD and 
naive fields and to [2,1] for a discussion of the 
constrained curve fitting we employ in extract- 
ing energies and matrix elements. A forthcom- 
ing paper [3] will give a detailed description of 
the methods used to obtain the results presented 
here. 



2. QUENCHED RESULTS 

The quenched configurations were generated 
using the tadpole-improved tree-level Symanzik 
action. Last year we reported difficulty in ob- 
taining reasonable fits for heavy-light correlators 
on a lattice with 1/a = 0.8 GeV when the 0{a?) 
staggered action (AsqTad) was used [1]. By the 
process of elimination, it was determined that 
the extended 3-link hopping of the Naik term in 
the temporal direction was responsible for the 
bad fits. This term is known to give rise to 
energetic negative-norm, ghost-like solutions to 
the free fermion dispersion relation. Since these 
states have energies proportional to l/fli, their 
effects on correlators should become negligible as 
at decreases. Indeed on an anisotropic lattice 
with 1/at — 3.7 GeV reasonable fits arc obtained 
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Table 1 



Summary of fits to pseudoscalar heavy-light correlators. 



volume action 


1/a, (GeV) 






EsUp = 0) (MeV) 


8^ X 201-link 


0.8 


1.0 


0.59 


735(10) 


8^ X 20AsqTad 


0.8 


1.0 


8.93 




8^ X 48AsqTad 


0.7 


5.3 


1.59 


790(36) 


8^ X 48AsqTad w/o 


f-Naik 0.7 


5.3 


1.03 


901(19) 


123 X 32 1-hnk 


1.0 


1.0 


0.48 


873(9) 


12^ X 32 AsqTad 


1.0 


1.0 


0.96 


765(9) 



with and without the temporal Naik term (see 
Table 1). Similarly, on an isotropic lattice with 
1/a = 1.0 GeV, there is no significant contamina- 
tion from negative norm states; a fit is shown in 
Fig. 1. 

We focus now on results on the 12^ x 32 lattice, 
where we used the AsqTad action with amo = 
0.10 and an ^(AQCD/Airo, a^)-improved NRQCD 
action with qMq — 5.0. The bare quark masses 
were estimated to coincide with the strange and 
bottom masses, respectively. Meson masses can 
be computed using either finite momentum cor- 
relators to construct the kinetic mass 

or using perturbation theory 

Mpert = E,i^{0) + ZMMo-Eo (2) 

where Zm is the heavy quark mass renormaliza- 
tion and £"0 is the self-energy constant. For this 
lattice and gMq = 5.0 we find ZmMq — Eq = 
Mo - 0.890 as + Mo(!?(a^). The result is Mpert = 
5.51 ± 0.45 GeV where the 0(a.'^^) uncertainty 
leads to the quoted error. The comparison of 
Mkin to Mpert IS showu in Fig. 2. 

Mass splittings can be computed directly from 
differences of simulation energies. We find Mb* — 
Mb^ = 25 ± 5 MeV, much smaller than the split- 
ting expected from the experimental Mb* —Mb = 
46 MeV, but in agreement with other quenched 
results. 

Matrix elements of the temporal QCD ax- 
ial vector current Aq are evaluated up through 
0{AQCD/Mo,as,as/aMo) by computing matrix 
elements of 



J, 



(0) 




Figure 1. Squares show the J^ '\t)J^"'{0) corre- 
lators on the 12^ x 32 lattice, and diamonds show 
the jQ°^\t)jQ^\o) correlator times —1. 



J. 



(1) 



2Mo 



975707 • 



(4) 



and using the operator matching relation 

Ao ^ (5) 

where the 0(Aqcd/Mo) contributions are ab- 
sorbed into 



T(l,sub) 
^0 



t'--^' ^, r 7'-"' 



(0) 



(6) 



With aMo = 5.0 we compute Za,, = 1 (0.208 ± 
0.003)as and Cio = -0.0997. Fits to appropriate 
local-local correlators yield 



(0|j(^'^"'^)|i?.) 



9 7570 Q 



(3) 



Applying (5) and 
(0|^o|S.) = SbMb, 



= - 0.034 ± 0.004 (stat). (7) 



(8) 
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Figure 2. Kinetic Bg mass on the quenched lat- 
tice. The sohd hues mark the range of Mpcrt and 
the dashed hne marks the experimental mass. 
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Figure 3. Kinetic Bg mass on the unquenched 
lattice. The dashed line marks the experimental 
mass. 



gives the quenched result 

fs, = 225 ± 9(stat) ± 20(p.t.) MeV. (9) 

The 20 MeV uncertainty is the estimate of the 
0{a1) error in Zaq- 

3. UNQUENCHED RESULTS 

Our work has just begun on the large suite of 
2+1 flavor MILC ensembles (see e.g. [4]). So far 
we have simulated the Bs system on the config- 
urations with a « 0.13 fm and dynamical bare 
masses auonii — 0.01 and auQiris ~ 0.05 (MILC 
absorbs a factor of uq into their definition of tuq 
compared to our convention). Our valence mass 
is auomo = 0.05. We have used the lattice scale 
determined from the T spectrum, 1/ar — 1.58(2) 
GeV [5]. This scale is a bit higher than the one 
quoted in [4], l/a^i = 1.46 GeV; however the 
heavy quark potential parameter ri is known only 
through quark models and so is not as rigorous of 
a quantity for setting the lattice spacing. 

In Fig. 3 we show the kinetic Bg mass, e.g. the 
resuh for p = (0,0, 1) gives Mkin(Bs) = 5.63(17) 
GeV. We find that the unquenching dramati- 
cally improves the B* — Bs splitting, bringing 
it up to 42.5 ± 3.7 MeV. The value we compute 
for 2Mb, - Mr = 1.361 (8) (17) GeV where the 
first error is dominated by the statistical error in 



Esira{Bs) and the second is due to the uncertainty 
in ax- Using ax implies that the strange sector 
actually corresponds to auouiQ w 0.04 [6]. From 
the physical Mb, — Mb we estimate our value 
of 2Mb, — M-f obtained with auQiriQ ~ 0.05 is 
too high by about 45 MeV, or 2Mb, - Mr = 
1.316(8)(17) GeV versus 1.278(5) GeV from ex- 
periment. Simulations with smaller auomo will 
permit a more rigorous interpolation to Mb, ■ 
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